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Abstract 

We study the three-body Coulomb problem in two dimensions and show how to calculate very 
accurately its quantum properties. The use of a convenient set of coordinates makes it possible 
to write the Schrodinger equation only using annihilation and creation operators of four harmonic 
oscillators, coupled by various terms of degree up to twelve. We analyse in details the discrete 
symmetry properties of the eigenstates. The energy levels and eigenstates of the two-dimensional 
helium atom are obtained numerically, by expanding the Schrodinger equation on a convenient 
basis set, that gives sparse banded matrices, and thus opens up the way to accurate and efficient 
calculations. We give some very accurate values of the energy levels of the first bound Rydberg 
series. Using the complex coordinate method, we are also able to calculate energies and widths of 
doubly excited states, i.e. resonances above the first ionization threshold. For the two-dimensional 
H~ ion, only one bound state is found. 
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I. INTRODUCTION 



Since the very beginning of quantum mechanics, the hehum atom has attracted much 
attention as it is one of the simplest system where the Schrodinger equation cannot be 
solved exactly. Recently, it has been understood that the lack of an exact solution is the 
direct quantum counterpart of the non-integrable character of the corresponding classical 
dynamics Indeed, it has been discovered that, for most initial conditions (positions 
and velocities of the two electrons), the classical dynamics is chaotic, with the total energy 
and the total angular momentum being the only constants of motion. Together with the 
development of sophisticated numerical methods for computing the quantum energy levels 
0, 1^, ^, 0], there have been major improvements on semiclassical techniques which 
allow to compute approximate values of the energy levels from the knowledge of the classical 
dynamics. The most dramatic success is the use of periodic orbit theory, where the energy 
levels are calculated from simple properties (action, period, stability...) of a (preferably 
large) set of classical periodic orbits 0] . Most of the quantum and semiclassical calculations 
concentrated on states with low total angular momentum for at least two reasons: firstly, 
these are the states experimentally prepared when using an optical excitation from a low 
excited state and, secondly, this is the situation where the classical dynamics is well known. 

Of special interest are the S states with zero total angular momentum. Classically, the 
motion of the two electrons takes place in a fixed plane. Thus, the classical dynamics is fully 
identical with the classical dynamics of the two-dimensional (2D) helium atom. It turns out 
that, although it seems to be a simpler system, there has been only very little interest in 
this 2D three body Coulomb problem and practically no quantum calculation. It is the aim 
of this paper to fill this hole. It can be also expected that, when a "real" 3D helium atom 
with low (or zero) initial momentum is exposed to an external perturbation, its response will 
not be very different from the one of the 2D atom, provided angular momentum does not 
play a crucial role in the physical processes involved. For example, when a helium atom is 
exposed to a strong non-resonant low-frequency electromagnetic field, it may absorb a large 
number of photons leading eventually to single or even double ionization. It seems likely 
that the correlation between the two electrons plays a major role in this process (especially 
in the generation of high harmonics of the electromagnetic field), while the total angular 
momentum remains relatively small. Another example is the production of doubly ionized 
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atoms where a process involving symmetric excitation of the two electrons (with zero total 
angular momentum) has been recently proposed |^. In these situations, the full 3D quantum 
calculation for such a system is not presently feasible, except for the very lowest states. On 
the other hand, a 2D quantum calculation seems reachable. This would allow to determine 
whether the proposed process is relevant or not. It is thus highly desirable to be able to 
compute accurately the quantum properties of the 2D helium atom. 

A second motivation to study the 2D three body Coulomb problem comes from semi- 
conductor physics. The study of excitons - the bound aggregate of an electron from the 
conduction band and a hole from the valence band, each particle with a given effective 
mass- is an important tool to study semi-conductors. In 1958, M.A. Lampert @ has shown 
that three body complexes called trions (an electron or a hole bound to an exciton) should 
be observable at low temperatures, and this was confirmed later by variational calculations, 
showing the stability of trions against dissociation into a exciton and a free electron or a 



hole (see [10| for references). Since, the progress in semiconductor technology have made 



possible the fabrication of quasi 2D systems. It was then realized [0, |Tl|] that in such 
systems, trions would have an increased stability due to the 2D confinement, and should 
thus be more easily observable. The trions are responsible for satellites on the excitonic lines 
in luminescence spectra. Several observations have been reported since the first one in 1993 
T^ , |TB|, |T^, |TB|, [TB|, |T^, and compared with theoretical predictions []TU|, |TB|. In this context. 



a precise calculation of the energy levels of the excitonic trions in a 2D system as a function 
of the ratio of the effective masses, with and without external field, is highly valuable, and 
justifies the methods and calculations introduced in this paper. The 2D hydrogen molecular 
ion H2 has also been studied in the frame of the Born Oppenheimer approximation in Ref. 
||T9| , where the first two electronic energy curves are given. 

The paper is organized as follows. In section II, we discuss the physical symmetries of 
the 2D three body Coulomb problem. We then introduce a new set of parabolic-like coordi- 
nates, give the expression of the Hamiltonian operator and show that we can find a basis in 
which the Schrodinger equation involves sparse banded matrices, allowing accurate numer- 
ical calculations. In section PTT] , we analyse the group structure of the discrete symmetries 
of the new Hamiltonian, showing that the complications introduced by the not one-to-one 
character of the change of coordinates can be taken into account exactly and actually does 
not lead to any difficulty. In section |V|, we first explain the detailed structure of the basis 
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Figure 1: The relative cartesian coordinates of particles 1 and 2 with respect to particle 3 are 
and (2:2, 2/2)- The interparticles distances ri, r2 and ri2. 

set that we use. We then discuss the structure of the expected energy spectrum in the case 
of a 2D helium atom with an infinite mass nucleus, and give the energies of the lowest levels 
in the bound Rydberg series, as well as - using the technique of complex coordinates - the 
energy and width of the first doubly excited resonance. 



II. THE SCHRODINGER EQUATION 



A. Hamiltonian 



The three body problem in two dimensions has 6 degrees of freedom that can be reduced 
to 4 in the center of mass frame. Here, as depicted in figure ^ ri and r2 denote the positions 
of particle 1 or 2 with respect to particle 3, and pi and p2 the conjugate momenta. In 
atomic units (such that h, 47reo, the mass m of the electron and the elementary charge are 
all equal to unity), the Hamiltonian writes, neglecting QED and relativistic effets: 

rr Pl^ , P2^ , P1-P2 , Q1Q3 , Q2Q3 . Q1Q2 /..x 

n = h 1 \ \ \ , (Ij 

2Aii3 2/i23 ri r2 ri2 

where is the mass of the third particle (in unit of the electron mass), and /ii3 (resp. /i23) 
is the reduced mass of particle 1 (resp. 2) and particle 3. Qi, Q2 and Q3 are the charges 
of the particules in unit of the elementary charge. ri2 is the distance between the particles 
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1 and 2. The 2D helium atom with a fixed nucleus corresponds to the case where is 
infinite, /^is = yU23 = 1, <5i = Q2 = -1 and Q3 = 2. 

As for the 3D three body problem we regularize the Schrodinger equation, i.e. remove 
the denominators, by multiplying it by 16 rir2ri2- The eigenstate |^) with energy E then 
satisfies the generalized linear eigenequation: 

16rir2ri2 + ^ + ) I*) + ^1^) = 16 nr2ri2 E |^), (2) 

\2/ii3 2/i23 TTL-i ) 

where 

= 16 (QiQs r2ri2 + QaQs riri2 + QiQs nra). (3) 

B. Symmetries 

The symmetries of the 2D three body problem are the rotational invariance around an 
axis (A) perpendicular to the plane, the parity 11 and, when particles 1 and 2 are identical, 
the exchange symmetry Pyi- In two dimensions, the parity operator 11 coincides with a 
rotation of angle tt around (A), so that 11 and the angular momentum are related by 

^=(-l)"^ (4) 

We also introduce the two commuting symmetries n^. (symmetry with respect to the x axis) 
and Ily (symmetry with respect to the y axis). They are related to total parity through 
Ti^y = Hyllx = n. The group generated by Ux, Hy and P12 is the so called D2h point 
group. It is an invariance group of the Hamiltonian (|l]), for identical particles 1 and 2. The 
symmetries Ux and II^^ both commute with parity, but not with the angular momentum 
since, for instance, 11^ = —LzUx- As a consequence, the eigenstates of the 2D three body 
Coulomb problem can be labelled by their angular momentum = 0, ±1, ±2 ... and by the 
exchange symmetry when particles 1 and 2 are identical. The spectrum corresponding to Ml 
and —Ml angular momenta are identical: this (Kramers) degeneracy is a direct consequence 
of the time reversal invariance of the problem Alternatively, the eigenstates could also 
be labelled by parity with respect to the x axis and the absolute value of the angular 
momentum. 

When the system is exposed to an external uniform electric field along the x axis, the 
angular momentum is no longer preserved. The only remaining symmetries are Ila; and P12. 
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C. Parabolic coordinates 



In order to perform efficient and accurate numerical calculations, we wish to obtain a 
sparse banded matrix representation of the linear problem (|^) where the non-zero matrix 
elements are known in a closed form. We thus have to find a basis set in which the various 
terms of the Hamiltonian have strong selection rules. This can be achieved for example if 
all terms of the Hamiltonian can be expanded in polynomial combinations of position and 
(conjugate) momentum coordinates: in such a case, the set of eigenstates of an harmonic 
oscillator is convenient. Our situation is slightly more complicated, because the Hamiltonian 
involves the interparticle distance. How to deal with such a problem is well known for the 
hydrogen atom: by introducing a set of so-called parabolic or semiparabolic coordinates 
||2^ , one can map the 2D hydrogen atom on an harmonic oscillator. The method used here 
for the 2D helium atom is inspired by such a treatment, although it is technically more 
complicated. 

If X and y are the cartesian coordinates of a point in a 2D space and z = x + iy is the 
associated complex number, the distance from the origin is r = \z\ = y/x'^ + y'^, and its 
expression involves a square root function. The square root can be removed if we introduce 
the complex variable Z = X + iY defined hj z = since r = = ^^^^^ . X and Y are 
the parabolic coordinates, related to x and y by: 

X = and y = XY. (5) 

The parabolic coordinates are extremely convenient to represent the hydrogen atom in two 
dimensions |^^, or the Stark effet of the 3D hydrogen atom [^. Of course, the correspon- 
dance between (X, Y) and (x, y) given in equation (|^) is not one to one. The difficulties 



related to that choice of coordinates are discussed in section HI B 



We now come to the case of three particles. The complex positions of particles 1 and 
2 with respect to particle 3 are zi and Z2, and Zi and Z2 are the associated parabolic 
coordinates. The interparticle distances then write ri = ^-^j-, r2 = ^-^j- and then ri2 = 
\zi — Z2\ = l(^i+^2)M(^i-^2)| ^ j£ introduce the two complex numbers Zp = and 
Zm = ^^^^ , the distance ri2 appears as the product of the moduli of Zp and Z^- Since we 
want to express ri2 using square moduli, we introduce a second parabolic transformation on 
both Zp and Z^ by setting Zp = ^ and Z^ = The three distances are then expressed 
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as the square moduli: 
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As a consequence, the three distances have polynomial expressions when they are expressed 
with the new coordinates {xp,yp,Xm,ym) defined by Ep = Xp + iyp and Em = Xm + iym- 
Those coordinates are related to the initial cartesian coordinates {xi,yi,X2,y2) by : 



Yg (-^p Vp "^^pUp "I" ym "^^mymj yp ~\~ "^Xpyp + X^ + IXmy-n 

— (^Xp 1/p + Xm Vrn) {Xpyp + Xmym) ; 
■^2 T77 {xp yp ~\~ '^Xpyp Xm ~\~ Vm '^Xmyrnj {xp 2Xpyp Xm + ~^ 2Xm,l/n 
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p 2/p -^m "I" 2/m,) {.XpUp XmV: 
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and the three distances are: 

'^i = ^ ({xp - Umf + iVp + Xmf) {{xp + ymf + (z/p - 

= ^ ((a;p + x^)^ + {yp + ((xp - Xmf + (z/p - 1/™)^) 
ri2 = j{xl + yf) (xl^ + yl 



(10) 



D. The Schrodinger equation 

The Schrodinger equation (0) can be written as: 
f T T T 1 



[ 2/ii3 2/i23 "^3 

where the kinetic energy terms are: 



^ ((Xp + Xmf + (Z/p + Z/m)^) ((Xp - X„)^ + (?/p - ymf) 



dxl dyl^ 

-\- ^[XpXm ypyrr 



Q2 \ 

- + + ■ dyl 



+ 



dxpdxm dypdym, 



2{xpy m ypXr, 



dxpdym dypdx 
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T2 = 

16 



K^m + Vm) yg^2 + 9y2 ) + l^P + ^/p^ [9^2^ + ^^2^ J 
- ZyXpXjYi ~\~ ypym) (00 ~^ ^) ^) ) ~^ ^y^pUm Vp^n 



ydxpdx„ oypdym) \dxpdym dypdxmj J 

^12 = f (a:' + ylY - {xl + yl,f) Uxl + yl) f + 1- ) - (x^ + yl) + ^ 



16 ^ ' V m ■ J yy ' ^ra, ^9^2 ' 9y2 j V P ' '^P' ^9^2^ 9y2^ 

'^i.XpXm ~l~ ypym) i^Xpy^n ypXm) \ {XpX^ -\- ypym) 



dypdxrn dxpdym, 



( d^ d^ 

\Xpym ypXm) 
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ydxpdx^ dypdymj j 

B = IQ riraria. (12) 

The expressions of B and V can be deduced from equations (|^) and (|TU|). The Jacobian 
of the coordinate transformation is 16 rir2ri2. The scalar product of two wave functions is 
given in appendix 0. 

The various terms in the Schrodinger equation ([TT|) are polynomials in the coordinates 
\Xpy ypj XjYi, y^ ) and their associated momenta (partial derivatives -id/d{xp,yp,Xra,ym})- 
The operators Ti, T2, T12, V and B can thus be expressed using the corresponding an- 
nihiliation and creation operators: 

= 71 + ^) ■ = 71 " ^) ■ 

This shows that the 2D three body Coulomb problem can be described using the annihilation 
and creation operators of 4 harmonic oscillators. The Hamiltonian is a polynomial of degree 
12 in the annihilation and creation operators. Consequently, it will be possible to choose a 
basis of tensorial products of Fock states of each harmonic oscillator, for which the operators 
involved in the Schrodinger equation exhibit strong coupling rules. 

From the annihilation and creation operators associated with the new coordinates, we 
introduce the right and left circular operators in the planes {xp, yp) and (x^, ym) defined by: 



«3 = {axm -^«s/m)/v^, (14) 
8 



Using the symbolic calculation language Maple V, we have calculated the normal ordered 
expression of the various operators involved in the Hamiltonian. Those expressions are too 
long to be published here. Indeed, the operators Ti and T2 contain 625 terms, T12 331, 
the potential operators riri2 and r2ri2 517, rir2 159 and B 1463. When particles 1 and 
2 are identical, the Hamiltonian involves the kinetic term Ti + T2 and the potential term 
(ri + r2)ri2 that have only 335 and 275 terms, because the terms of Ti and T2 that do not 
commute with the exchange operator P12 cancel out. 



E. Angular momentum 

The angular momentum Lz has a very simple expression when expressed with the 
{xp,yp,Xm,ym) coordinates: 

( d d d d \ 

Lz = -i\ xi- yi- h X27^ 1/2- 



dyi dxi dy2 dx 



2 , 



T - ^ ( ^ ^ ^ ^ \ 

4 \ dyp '^dxp dym oxm) 

The relation z = ^ between the cartesian and the parabolic complex numbers shows that a 
rotation of angle ^ on Z is a rotation of 29 on z. Consequently, a factor of 2 appears in the 
expression of the angular momentum in parabolic coordinates p2|. Since we have performed 
two successive parabolic transformations to define the (xp, ?/p, Xm, y-m) coordinates, we have a 
factor 4 in the denominator of equation (p^. With the annihilation and creation operators 
(p!^), the angular momentum simply writes : 



Lz = {N^-N2 + N;-N,)/A, (16) 

where the number operators are Ni = afoj. They are related to the number operators 
corresponding the the annihilation and creation operators , ... given in eq. (|13D by: 



N, + N, = N,^ + Ny^. (17) 
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III. DISCRETE SYMMETRIES 



Xp 


Xp 


ITy . Xp 


Vp - 


-Vp 


Vp 


Xm 


Xm 


Xm 


ym ~ 


Uin 


ym 



A. Physical symmetries 

The Hamiltonian (0) has two discrete symmetries, n^, and IIj^, which are the symmetries 
with respect to two orthogonal axis in the physical plane. Using the new (xp, ?/p, Xm-, ym) 
coordinates, they can be expressed for instance as: 

• Xp > Xp iiy . Xp > (-^p ~i~ yp) I 

[xp-yp)l^pl ^^^^ 
[xm + y m 
' {xm - ym)/y/^- 

Moreover, if particles 1 and 2 are identical, the Hamiltonian commutes with the exchange 
operator P12. The effect of P12 on the (xp, x^, ym) coordinates is: 



(19) 



Obviously, the Schrodinger equation ( ]Tl| ) written with the (xp, j/p, x^, ym) coordinates is 
invariant under these transformations. 



P12 '■ Xp — 
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B. "Additionnal" symmetries 

In this section, we analyse the constraints that the physical wave functions must satisfy. 
We first recall what happens in the case of a single parabolic transformation. The parabolic 
transformation {X, Y) —>■ (x, y) defined in equation is a one-to-one mapping of the 
quarter of plane {X > 0,Y > 0) onto the half-plane {x,y > 0). Here, the transformation 
is used to represent the full cartesian plane (x, y) by extending the domains of X and Y 
to ] — 00, +oo[. That way, we obtain a double mapping of the cartesian plane since {X, Y) 
and {—X, —Y) are mapped on the same point. Consequently, the Hamiltonian written with 
the parabolic coordinates has a new discrete symmetry {X, Y) —>■ {—X, —Y), i.e. the parity 
with respect to (X, Y). The physical wave function must be a single- valued function of the 
initial coordinates {x,y) i.e. must fulfill \E'(X, F) = \E'(— X, —Y). Any function of \E'(X, y) 
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which satisfies the Schrodinger equation written in the (X, Y) coordinates but does not obey 
the constraint Y) = X, —Y) is to be rejected as an unphysical solution. 

In the particular case where the wave function is expanded on a basis built with tensorial 
products of harmonic oscillator eigenstates: 

1^) = E Cn^,ny \nx) ® (20) 
nx,n-Y 

the physical wave function expansion of equation (pO]) is restricted to the even values of 



nx + ny, because the parity of the Fock state \n) is (—1)" |22 . 

This property can be extended to the case of the transformation given in equation 
that give the cartesian coordinates versus the new coordinates {xp,yp,Xm,ym)- Because we 
perform four parabohc transformations to obtain the {xp,yp,Xm,ym) coordinates from the 
initial cartesian coordinates, there are four "additionnal" discrete symmetries which leave the 
Schrodinger equation ([TI|) invariant. We denote them Hi defined as (Xi, Yi) (— Xi, — Yi), 
112, Hp and n^^. The effects of those symmetries on the {xp,yp,Xm,ym) coordinates are: 

n^i • Xp > y^Yi II2 . Xp > Xfyi 
yp Xm yp ^ ym 

m yp •'^m -^p 

Vin ^ Xp ym ^ Vp 



(21) 



rip . Xp ^ Xp rirn • ^ 

yp ^ ~yp yp ^ yp 

Um ^ ym Um ^ Um- 



C. Symmetries of the wave function 

The group G generated by the lix, Hj,, P12 and the Hi, 112, Hp, 11^ symmetries is an 
invariance group of the Schrodinger equation (0). It is studied in details and its character 
table is given in appendix 0. 

In order to be singlevalued in the geometrical space (xi, yi, X2, 2/2), the wave function 
\E'(xp, Xm, must be invariant under any "additional" symmetry introduced by the 
non one-to-one change of coordinates, i.e. under any of the transformations Hi, 112, Hp, 
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Um- Then, the wave function must belong to an irreductible representation of G for which 
the character of any "additional" symmetry is equal to its dimension. There are only 8 
representations with this property, all being one- dimensional, that correspond to the first 8 
lines of the character table given in appendix Consequently, the physical eigenfunctions 
\l/(xp. Hp, Xm, ym) can be distinguished only by their symmetry properties with respect to Ha,, 
Ily and Pi2. The 8 physical irreductible representations of G are those of the group D2h (or 
mmm), of order 8, already mentionned in section ( |1I B]) . The application that maps each 
"additional" symmetry on the identity is a group homomorphic mapping of G on D2h- 

Finally, we have shown here that all energy levels belong to a one-dimensional represen- 
tation of the discrete symmetry group of the Schrodinger equation, and are thus expected 
to be non degenerate (except for the {M^, —Ml) mentionned above). Moreover, using the 
(xp. Up, Xm, Um) coordinates does not introduce extra representations which cannot be distin- 
guished from the physical ones. The wave functions can be described using a basis exhibiting 
the relevant symmetry properties with respect to n^,., 11^ and P12, or and P12. The second 
feature will be extensively used in the numerical implementation. 

In other words, among all solutions of the Schrodinger equation (pJ^ in the (xp, yp, Xm, ym) 
coordinates, sorting out the unphysical solutions is rather straightforward and one is left only 
with the physical symmetries of the initial system. 



IV. NUMERICAL SOLUTION 



A. Basis set 



1. Basis structure 



To perform numerical calculations of the eigenenergies and eigenstates of the three body 
Coulomb problem, we expand the Schrodinger equation on a convenient basis, and then solve 
a linear eigenvalue problem. Because the different terms of the Hamiltonian have polynomial 
expressions in the annihilation and creation operators, we obtain strong selection rules if 
we choose basis functions that are tensorial products of Fock states \ni) of the harmonic 
oscillator described by the circular annihilation operator Oj. Namely, we set : 

\ni,n2,n3,n4,) = \ni) ® \n2) ® \n-i) ® 1^4). (22) 
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The indices are then positive integers. The basis functions are eigenf unctions of the 
angular momentum, corresponding to the integer eigenvalue: 

Ml = (ni-n2 + n3-n4)/4. (23) 

The wavefunctions of these basis states are simple. Indeed, they are just eigenstates of an 
harmonic oscillator along the various coordinates. In the {xp, yp, Xm, ym) coordinates, they 
should appear as products of Hermite polynomials and Gaussian functions of the coordinates. 
As we use circular creation-annihilation operators, equation ([T^), the associated eigenstates 
of the two-dimensional harmonic oscillators in the {xp,yp) and {xm,ym) planes are easily 
written in polar coordinates as the product of a exp{i(j)) term with an exponential and a 
Laguerre polynomial of the squared radius. The explicit expressions of such states can be 
found in [^]. 

We have previously shown that the two successive parabolic transformations introduce 
"additional" unphysical states. The physical solutions can be selected using a basis set that 
is even with respect to all the "additional" symmetries. This choice is performed in two steps. 
First, both ni+n2 and 71,3+714 have to be even numbers. Indeed, from equation (|l3), ni+n2 = 
^xp + ny^ and ^3 + ^4 = n^^+Uy^ and the even representations for Up and Um corresponds 
to even values of nx^+riy^^ and nx^+riy^. Secondly, because the transformation (1, 2, 3, 4) — »• 
(3,4,1,2) on the annihiltion and creation operators commutes with the Hamiltonian and 
corresponds to the identity in the physical space, the basis functions have to be chosen as 
the symmetric combinations: 

|ni, n2, ns, = |ni, n2, n^, n^) + |n3, 114, ni, n2). (24) 

Of course, this symmetrised state remains an eigenstate of the angular momentum, with the 
same eigenvalue Ml- Taking into account the even parity of ni + n2 and n^ + n^, and thus of 
rii — n2 and —n^, and the expression of Ml, we obtain that ni — n2 (mod 4) and 123 — 
(mod 4) are simultaneously equal to either or 2. We then set: 

C12 = (ni - n2) (mod 4) = (723 - ua) (mod 4). (25) 

When particles 1 and 2 are identical, the Hilbert space can be split into a singlet sub- 
space corresponding to C12 = 0, and a triplet subspace corresponding to Cu = 2. Here, 
singlet means symmetric with respect to the exchange operator P12 whereas triplet means 
antisymmetric. 
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We can now define precisely the basis set corresponding to the physical states with angular 
momentum Ml and either singlet or triplet exchange symmetry. Since the quadruplet of 
indices (ni, n2, ns, 714) and (ns, 714, ni, give the same symmetrised ket in equation (p^D, 
we have only to consider one of the two quadruplet to label uniquely the symmetrised basis. 
Consequently, for singlet states, we set: 

C12 = 0, (^1 > or {rii = and n2 > ^4)} , (26) 

and for triplet states: 

^anti-sym ^ ||^^^ ^ ^ - n2 + Us - = AMl, Hi > 0, 

Ci2 = 2, (ni > ns or (ni = and n2 > rii)} . (27) 



2. Selection rules and matrix elements 



Two basis vectors \ni, n2, n^, and \ni + 5ni, n2 + 5?7,2, + 5n^, 71,4 + 577,4) are coupled 
by the Hamiltonian if the shifts 5ni correspond to one of the 225 allowed coupling rules. 
Because the Hamiltonian commutes with the total angular momentum, they all obey 5ni — 
5n2 + — drii = 0. Among them, 159 rules preserve the exchange symetry while 66 do 
not. The 159 rules that appear for the operators Ti +T2, T12, (ri + r2)ri2, rir2 and B, obey 
6ni — 6n2 = —{Sn^ — Sn^) = or 6ni — 6n2 = —{Sn^ — Sn^) = ±4, and are shown in figure |^. 
The 66 ones verify 6ni — 6n2 = —{Sn^ — 6n4) = ±2. They appear if the exchange symmetry 
is broken (mi 7^ m2 or Qi 7^ Q2) in the kinetic terms Ti, T2, and the potential terms riri2 
and r2ri2. 

Since the Hamiltonian has been written in normal order, the derivation of the matrix 



elements is straigtforward. They are too numerous to be written explicitely here ||30|. We 
only give two matrix elements of the kinetic operator of the 2D helium Ti + T2 between two 
unsymmetrised basis vectors: 

(ni,n2,n3,n4|(Ti + T2)|ni, 772,773,^4) = (28) 

1/4 (721 + n2 + 1) (7Z3 + 7Z4 + 1) 

(ni^ + 4 ?7i 7l2 - Til 773 + ni ?74 + 7l2^ + 712 173 - '"-2 
+?73^ + 4 ?73 7Z4 + ?74^ + 3 Til + 3 722 + 3 773 + 3 774 + 

14 



6n. 



(§)(§)( • )(§) 

m (§) ® ® ® 

(§) ® o o o 

(§) o o o o 

(§) o o o o 

($) o o o o 



6n, 



• ® 

(§) 

(§) 
® 



Figure 2: The 159 selection rules that preserve the exchange symmetry are depicted in the (ni, 713) 
space. The dark circles corresponds to the 61 rules 5ni = 6n2 and ^^3 = 6n4, the full line circles 
to the 49 rules 5ni — 5n2 = 4 and Sn^ — dn^ = —4 and the dashed line circles to the 49 rules 
6ni — 6n2 = —4 and dn^ — Sn^ = 4. 

as well as the matrix element corresponding to the selection rule 5i = ^2 = ^3 = ^4 = 1: 

(ni + l,n2 + l,n3 + l,n4 + l|(Ti + T2)|ni, n2, 7^3, ^4) = (29) 



-1/4 V^4 + Wns + Wn2 + Wni + 1 

(jii^ + 5 rii n2 — 2 rii + n2^ — 2 r;,2 ^3 + n^^ + 5 713 714 

+71^ + 5 77-1 + 5 77-2 + 5 77-3 + 5 n4 + 12^ . 



3. Numerical implementation 

For the numerical calculations, we have chosen to truncate the basis defined by equation 
(PB|) or ( P7D using the condition ni + n2 + 773 + 774 < Nhase- Because the angular momentum is 
fixed, we have only 3 independant indices, and the size of the basis is roughly iV^^^g/192. The 
basis B is then ordered in order to represent the Schrodinger equation using band matrices 
as narrow as possible. The eigenvalue problem is then solved using the Lanczos algorithm 
which makes it possible to compute a few eigenvalues in the range of interest. 
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4- Variational parameter 



So far, the natural length scale of the problem is the Bohr radius oq. Because it is not 
necessarilly the best suited one, we introduce the length scale as a~^ao. The Schrodinger 
equation ([TT|) writes: 

When the basis is truncated, the length scale a becomes a variational parameter (i.e. the 
calculated energy levels should not depend on a is the basis set is large enough) that has 
to be numerically optimized. All the numerical results presented in this paper are obtained 
with a close to 0.4. All the digits of the energy levels given in the Tables are significant. 
The uncertainty on the results is thus 1 on the last figure, and the relative accuracy reaches 
the 10-^3 level. 



B. The 2D helium atom without electron interaction 

Let us consider the 2D helium with a fixed nucleus of charge Q3 = 2 (the mass is 
infinite). The Schrodinger equation equation ( pOD simply writes : 

ja^^llZl + ^syj 1^) = a'^ EB\^), (31) 

where V = —32 (ri + r2)ri2 + 16rir2. If the 16rir2 term in the potential energy is removed, 
the three body problem corresponds to two independent 2D hydrogen atoms with a nucleus 
of charge Q = 2. The spectrum of the 2D hydrogen atom is well known, and is given by the 
series fB^: 



2{N -l/2f 

where > 1 is the principal quantum number and — A^ + 1<M<A^ — 1 the angular 
momentum of the electron; the degeneracy is 2A^— 1. The structure of the energy spectrum is 
very similar to the 3D energy spectrum, the only difference being that the effective quantum 
number A^ — 1/2 is a half-integer ranging from 1/2 to infinity rather than a non-negative 
integer. 

For the helium atom without electronic interaction, the spectrum is thus given by: 

4 4 

En.,m. + E^.^M, = -2(Ar,_ 1/2)2 ' 2{N2-l/2r ^^^^ 
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where A^i and N2 are the principal quantum numbers of the two electrons. The essential 
degeneracy is 2(2A^i-l)(2A^2-l) if A^i 7^ N2 and (2iVi-l)2 otherwise 0. The total angular 
momentum is simply given by = Mi + M2. The states of total angular momentum 
correspond to the indices {Ni, Mi, N2, M2) and {N2, M2, Ni, Mi). Those degenerate states 
give symmetric (singulet) and antisymmetric (triplet) states when the two quadruplet are 
different and only one symmetric state if they are equal. Finally, the energy levels can be 
labelled by Ni,N2, Ml and Pu- The degeneracy of this configuration is given by the number 
of solutions of Ml = Mi + M2 taking into account the boundaries on Mi and M2. 

By solving the Schrodinger equation for an angular momentum between -3 and 3, and 
for the two exchange symmetries, we have checked that our method gives the expected 
eigenenergies and degeneracies. 

We have then checked the effect of the electronic interaction by introducing it perturba- 
tively as e/ri2. We have numerically computed the ground state energy of the three body 
problem as a function of e and observed a linear behaviour, as expected from first order 
perturbation theory. The slope in atomic units is 4.70(1), in agreement with the slope 37r/2 
predicted by first order perturbation theory (see appendix 0) . 

C. The 2D helium atom 

The l/ri2 term describing the electronic repulsion is now taken into account. This does 
not affect the positions of the various ionization thresholds (as the electron interaction 
vanishes at large distance). There is an infinite number of single ionization thresholds 
associated with the principal quantum number of the hydrogenic state of the resulting He"*" 
ion, given by energies: 

4 

= "2(iV- 1/2)2 • (34) 
These single ionization thresholds form a series which converge to the double ionization 
threshold at zero energy. 

Consequently, one expects bound states below energy Ii = —8 a.u., resonance (doubly 
excited states) between Ii and zero, and only continua above. 
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1. Bound states 



The lowest energy levels of the 2D helium below the first ionization limit are given in 
Table | for the singlet states and in Table ^ for the triplet states. For each value of M^, we 
obtain a Rydberg series converging to the = 1 threshold. For such excited states, the outer 
electron lies far from the nucleus while the inner electron is essentially in its ground state and 
lies very close to the nucleus. Because this picture gives two very different roles to the two 
electrons, it results in a new set of quantum numbers, namely (A^, M) for the inner electron 
and (n, m) for the outer one. A brutal but useful approximation is to neglect the effect of 
the outer electron on the inner one, i.e. consider that the inner electron in the hydrogenic 
state A^ = 1, M = while the outer electron sees a point charge Q = 1 (the charge +2 of 
the nucleus screened by the charge -1 of the inner electron) at the origin, resulting in an 
energy spectrum —8 — l/(2(n — 1/2)^), where n is the principal quantum number of the 
(hydrogenic) outer electron. This is of course only an approximation. Deviations from it 
can be measured through the quantum defect 5n,m defined directly from the energy levels 
through: 

= " 2(n - 1/2 - V™) • ('^^ 

If the previous approximation were exact, the quantum defects will all be zero. Hence, 
deviations from zero and evolutions with n and m directly measure the breaking of the 
approximation. The results shown in Tables | and |I| show that - alike the 3D helium atom 
- the quantum defect in a given series tend to a constant value as n — > oo. When |m| is 
increased, the outer electron is repelled from the nucleus by the centrifugal energy barrier 
and fills less the presence of the inner electron. It is thus expected that the quantum defects 
will decrease with increasing \m\ and this is fully confirmed by our "exact" diagonalizations, 
see Tables |, ^ and |ITT[ Also, in the triplet states, the wave function in configuration space 
is antisymmetric, so that the two electrons cannot seat at the same place. Consequently, 
they feel each other less efficiently, resulting in a lower interaction energy and consequently 
a larger quantum defect. Again, our "exact" calculations confirm this behaviour. 

Finally, it is interesting to compare the 2D and 3D situations. Although the Rydberg 
series are similar in both cases, this is not true for the ground state. Indeed, the binding 



energy of the inner electron in its ground state is 8 a.w., see equation (|3^) , in 2D, that is 
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A^, M, n, m Energy (a.u.) (5, 



1,0, 


1,0 





-11.899 822 342 953 


0.1419 


1,0, 


2,0 





-8.250 463 875 379 


0.0871 


1,0, 


3,0 





-8.085 842 792 777 


0.0866 


1,0, 


4,0 





-8.042 911 Oil 139 


0.0865 


1,0, 


5,0 





-8.025 668 309 76 


0.0864 


1,0, 


6,0 





-8.017 061 08 


0.0864 


1,0, 


2,1 


1 


-8.211 542 089 886 


-0.0374 


1,0, 


3,1 


1 


-8.077 637 328 985 


-0.0378 


1,0, 


4,1 


1 


-8.039 947 878 


-0.0378 


1,0, 


5,1 


1 


-8.024 280 94 


-0.0379 


1,0, 


3,2 


2 


-8.079 805 619 119 


-0.0030 


1,0, 


4,2 


2 


-8.040 745 817 


-0.0030 


1,0, 


5,2 


2 


-8.024 657 76 


-0.0031 


1,0, 


6,2 


2 


-8.016 51 


-0.0031 


1,0, 


7,2 


2 


-8.011 80 





Table I: Energy levels of the singlet states of the 2D helium atom (with infinite mass of the nucleus), 
below the first ionization threshold. The optimum variational parameter a is close to 0.4. For most 
of the states, the basis truncation is given by N^ase = 200. The basis size is then 43626 for singlet 
Ml = states, and slightly decreases with Ml. For the (1, 0, 4, 0), (1, 0, 5, 0) and (1, 0, 6, 0) we use 
^base = 240 and a basis size of 74801. In the fourth column, 5n,m is the quantum defect of the 
state, as deduced from equation (|35|). 

four times more than in the 3D helium atom. Almost the same ratio 4 is observed between 
the total binding energy of the ground state: 11.90 a.u. in 2D versus 2.91 a.u. in 3D. 

For singly excited states, the core is also 4 times smaller in 2D. We then expect a smaller 
core penetration because the centrifugal barrier is almost the same in the 2D and 3D systems, 
and also a smaller core polarisation by the outer electron, resulting in smaller quantum 
defects in the 2D case. The comparison of the 2D and the 3D quantum defects in Table |T| 
is consistent with this interpretation. 
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A^, M, n, m 


Ml Energy (a.u.) 




1,0,2,0 





-8.295 963 728 090 


0.2002 


1,0,3,0 





-8.094 583 618 582 


0.2008 


1,0,4,0 





-8.045 941 305 572 


0.2010 


1,0,5,0 





-8.027 055 169 


0.2011 


1,0,6,0 





-8.017 807 


0.2011 


1,0,2,1 


1 


-8.225 772 173 259 


0.0118 


1,0,3,1 


1 


-8.080 919 691 737 


0.0142 


1,0,4,1 


1 


-8.041 165 882 92 


0.0149 


1,0,5,1 


1 


-8.024 858 500 


0.0152 


1,0,3,2 


2 


-8.079 819 688 304 


-0.0028 


1,0,4,2 


2 


-8.040 751 693 48 


-0.0028 


1,0,5,2 


2 


-8.024 661 158 


-0.0028 


1,0,6,2 


2 


-8.016 512 


-0.0028 



Table II: Energies of the triplet states of the 2D helium atom (with infinite mass of the nucleus), 
below the first ionization threshold. The optimum variational parameter a is close to 0.4. The 
basis truncation is given by Nf,ase = 200. The basis size is 43550 for triplet Ml = states. In the 
fourth column, 5n,m is the quantum defect of the state, as deduced from equation (|35|). 



3D Rydberg series 


2D Rydb 


erg series 




6 


N,Ml 


6 




0.140 


1,0 


0.0864 


Sge 


0.299 


1,0 


0.2011 


Ipo 


-0.012 


1,1 


-0.0379 


3po 


0.068 


1,1 


0.0152 


lp)e 


0.0021 


1,2 


-0.0031 


3]je 


0.0028 


1,2 


-0.0028 



Table III: Quantum defects for various series of the 2D and 3D helium atoms below the first 
ionization threshold. The values in the 3D case are calculated from the energies given in p^]. The 
values in the 2D case are the limits of 6n,m for large values of n. 
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2. Resonances 



Above the first ionization threshold, the spectrum contains resonances embedded into the 



continuum. They can be numerically separated using the complex rotation method |2^, |26 
(also known as the method of complex coordinates), where the positions and momenta 
r and p are reespectively changed into re*^ and pe^*^. Here, it is simply implemented 
using a complex length scale a = |a|e*^ (in order to preserve their canonical commutation 
relations) This results in a "complex rotated" non hermitian Hamiltonian whose eigenvalues 
are complex. In the complex energy plane, the resonances do not depend on the angle 9 
while the continua are rotated by an angle of 26 around the ionization thresholds. The first 
resonance of the 2D helium atom (infinite mass of the nucleus) is obtained for zero angular 
momentum and singlet exchange symmetry. Its energy is: 

E = -1.411 496 328(1) - i 0.001 241 734(l)a.u. (36) 

It is obtained for a rotation angle 9 ^ 0.4, a length scale a ~ 0.35, Nhase = ^50, and a basis 
size of 18696. 

The energy structure of the resonances is illustrated in the case of the singlet Ml = 
states in figure ^. The electronic repulsion is included in the potential energy as e/ri2 
with < e < 1. One can follow the energy levels as a function of e from the independent 
electron case (e = 0) to the helium case (e = 1). For e = 0, the levels correspond to 
the Rydberg series N = 2, n > 2 converging to I2 = —8/9 a.u.. The degeneracy of the 
(A^ = 2,n = 2) configuration is 9, with 3 states of zero total angular momentum. Two of 
them have the singlet symmetry and one the triplet symmetry. For n > 2, the degeneracy 
of the configuration is 18, with 6 states corresponding to Mi = (3 singlet and 3 triplet 
states). Consequently, for e = 0, the first Ml = singlet resonance is doubly degenerate, 
and the following ones are triply degenerate. The introduction of the electronic interaction 
removes the degeneracy. 

D. The 2D H~ ion 

The H~ ion with a fixed nucleus is obtained by setting = 1. We obtain only one 
bound state below the first ionisation limit (at -2 atomic units from equation (P^)), with 
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Figure 3: Energy levels of the singlet Ml = resonances of the 2D helium atom (nucleus with 
infinite mass) between the first (-8 a.u.) and the second (-0.888. . . a.u.) ionisation thresholds as a 
function of the magnitude e of the electronic interaction included in the potential energy as e/ri2. 
The right part of the figure is a zoom close to the second ionisation threshold. All the energies 
presented on this picture are well converged and obtained with a = 0.5, 6 = 0.12, Nf,ase = 150. 
The basis size is 18696. Because the method of complex coordinates is used, these energy levels 
are not bound states, but resonances. At the scale of this figure, their widths are very small. 

zero angular momentum and singlet exchange symmetry. Its energy is: 



It is obtained for a ft; 0.4, Nhase = 220, the basis size being 57820. 
V. CONCLUSION 

We have introduced a new set of coordinates to represent the 2D three body Coulomb 
problem and given the resulting Schrodinger equation. We have discussed the discrete 
symmetry group properties of the equation and shown that it can be numerically solved 
very efficiently, using a convenient basis set for which the Schrodinger equation involves 
sparse banded matrices. The convergence of the calculations is very good and the numerical 
results are extremely accurate. This is demonstrated in the case of the 2D Hehum atom 
(with inifinite mass of the nucleus) , for which the lowest energy levels in the bound Rydberg 
series are given with a relative accuracy in the 10"^ to 10~^^ range. 



E = -2.240 275 363 589(l)a.u. 



(37) 
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The method developped in this paper provides an efficient tool for studying the dynamics 
of the 2D Hehum atom in an external electric field aligned along the x axis. Indeed, with a 
field of strength F, we must add the external potential energy term Ve^t = 16 rir2ri2{xi + 
X2) F to equation (^. The only remaining symmetries are then the exchange symmetry P12 
and the symmetry with respect to the x cartesian axis Tlx. In such a case, the convenient 
basis set can be defined, from equation (P^D, by: 

1^1, n2, ris, n^Y = \ni, n2, n^, n4)+ + e \n2, rii, n^, ^3) + , (38) 

with e = ±1 for even or odd states with respect to Ila;. Because V^xt is a polynomial in 
the new coordinates, it exhibits selection rules, making an accurate diagonalisation of the 
Schrodinger equation still possible. 

The motion of the nucleus can easily be taken into account, including the T12 contribu- 
tion to the Hamiltonian in equation ([TT|) . That way, it will be possible to determine very 
accurately the ground state energy of excitonic trions, as a function of the electron to hole 
mass ratio. This work is in progress. 
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Appendix A: CHARACTER TABLE 

The discrete symmetry group G of the Schrodinger equation (|rT|) written in the 
{xp,yp,Xm,ym) coordinates is studied here. It is generated by the 11^,., II^^, P12 physical 
symmetries and the Hi, 112, Hp, 11^ "additional" symmetries, defined in equations (|T8|), 
(p!9|) and (pi]) . Its structure has been studied following standard methods of the theory of fi- 
nite groups (see for example ^). Because all the generators of the group can be seen as 



a permutation among the sixteen quantitites Xp, ?/p, x^, {xp + Vp) / (xp — yp) /\/2, {xm + 
ym)/y/^, {xm — ym)/V^ and the opposite values, the G group appears as a subgroup of the 
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permutation group of sixteen elements. Thanks to this property, it has been easily studied 
using the permutation group package provided by the Maple language. 

The group G contains 128 elements, in 29 classes. It has 16 one-dimensional, 8 two- 
dimensional and 5 four- dimensional irreductible representations. Table represents its 
complete character table. It has been obtained using the method described in [Q. Let the 
classes of G be Ki with 1 < i < n. The set KiKj, set of the products of any element of the 
class Ki by any element of the class Kj, is made of complete classes. Calling Ciji the number 
of occurrences of the class Ki in the product KiKj, one can symbolically write: 

K,Kj=J2ci,iKi. (Al) 
I 

This property is used to obtain relations between characters xl of the class Ki in an 
irreductible representation R : 

29 

gigj Xi Xj - Xe 2^ (^iji 9iXi ■ [■^^i 
1=1 

where , Qi is the number of elements of the class Ki and Xe is the character of the identity E, 
i. e. the dimension of the representation R. The n characters of an irreductible representation 
appear to be the solutions of the n(n+ l)/2 quadratic equations obtained from ([A2| ) for any 
couple of (i, j). Then, to construct the character table of the group, three steps are necessary: 
first, the group has to be separated in classes, and the number Qi are obtained. Second, the 
numbers Ciji are computed, and last the system of equations ^ is solved. Obviously, there 
n different sets of solutions, corresponding to the n irreductible representations. 



Appendix B: 

The scalar product of two wave functions |\E'(^)) and |\E''^^)) in the {xp,yp, Xm,ym) coordi- 
nates is : 

^^(i)|^(2)^ _ _L JJJJ xi,(^)[xp,yp,Xm,ym)* B ^^'^\xp, y^, x^, y^) dxp dyp dxm dy^, (Bl) 

where B is given by equation (|T^). The integrals are calculated from — oo to +oo. The 
factor 1/16 = 1/2'^ comes from the four double mappings of the space introduced by the 
change to parabolic coordinates. 
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We now calculate the average value of l/ri2 for the ground state |^I/o,o) of a 2D helium 
atom without electronic repulsion, that is: 

1 



a 



0,0 



ri2 



0,0 



(B2) 



The normalised wave function of the ground state of a 2D hydrogenic atom with a nucleus 
of charge Q is : 



so that: 



^o,o(n,r2) = ^^0,0(^1) ^0,0(^2) 



8Q^ 



-2Q{ri+r2) 



TT 



(B4) 



We now evaluate a using the (xp, x^, ym) coordinates. Since the jacobian of the coordinate 
transformation is -B = 16rir2ri2, a writes: 



a 



16rir2|\l'o 



(B5) 



where ri and r2 are given by equations (|ToD. The integrals are calculated from —00 to +00. 

To evaluate a, we represent the {xp,yp) and {xm,ym) planes using polar coordinates 
(rp,9p) and {rm,Om) and obtain from equation ([T0|): 



n = ^{r'p + ri + 2rlrlcos{2ep-2e„ 
r2 = ^ « + ri- 2rlrlcos{2ep - 29„ 



ri2 



16 

p m 



(B6) 



The ground state wave function is then: 



-Q{4+rt)/2 



TT 



(B7) 



so that : 



(rj + + 2ry^ (l - 2005^(2^^^ - 26^))) e-^^^v^^'rn)!^ ^^^^^ ^^^^^ 



(B8) 

The integration over Qp and Q^a gives for the angular dependant term and 47r^ for the 
independant one. The integration over Vp and involves Gaussian integrals that give 

Svrg 



(B9) 
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that is 37r/2 when Q = 2. In the 3D case, a is evaluated to 5Q/8 in p9[. The ratio 
o"2d/o"3D = 67r/5 ~ 3.77 is close to 4, because the 2D ground state wave function is four 
times smaller than the 3D wave function. 
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Table IV: Character table of the discrete symmetry group G of the Schrodinger equation (11). The classes and the irreductible representations 
have been organized in order to obtain the character table of the group in the upper left corner (in bold figures). The last three classes 
are those of the "additional" symmetries 112, Hi and Hp {Hp and IIj^ belong to the same class). The first line gives the number of elements 
in each class. We finally give one element of each of the 29 classes : E, Hy, H^, Hyllx, P12, HyPu, U^Pu, HyH^Pu, 11^11^, n2nj;n^, n2Pi2, 

U1U2, n2n^Pi2, UpUyU^Pu, UpPi2, U2U^, U2Uy, U2UyU^, niHan^, n2n2,n^, U2UyU^Pu, ii2UyPi2, n^^n^n^, UpUyU^Um, UiUyU^p^, 

IlyIlxPi2^, II2, Hi, Hp, some of them are reported in the second line of the Table. 



